Calculations of polarizabilities and hyperpolarizabilities for the Be+ ion 

Li-Yan Tang^'^, Jun-Yi Zhang^, Zong-Chao Yan^'^, Ting-Yun Shi^, James F. Babb^ and J. Mitroy^ 
^ State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, 
Chinese Academy of Sciences, Wuhan 430071, P. R. China 
'^Graduate School of the Chinese Academy of Sciences, Beijing 100049, P. R. China 
'^School of Engineering, Charles Darwin University, Darwin NT 0909, Australia 
^ Center for Cold Atom Physics, Chinese Academy of Sciences, Wuhan 430071, P. R. China 
^Department of Physics, Wuhan University, Wuhan 430072, 
China and Department of Physics, University of New Brunswick, 
Fredericton 4400, New Brunswick, Canada and 
^ITAMP, Harvard- Smithsonian Center for Astrophysics, 
60 Garden St., Cambridge, MA 02138, USA 
(Dated: September 24, 2009) 

The polarizabilities and hyperpolarizabilities of the Be+ ion in the 2 ^S state and the 
2^P state are determined. Calculations are performed using two independent methods: i) 
variationally determined wave functions using Hylleraas basis set expansions and ii) single 
electron calculations utilizing a frozen-core Hamiltonian. The first few parameters in the 
long-range interaction potential between a Be+ ion and a H, He, or Li atom, and the lead- 
ing parameters of the effective potential for the high-L Rydberg states of beryllium were 
also computed. All the values reported are the results of calculations close to convergence. 
Comparisons are made with published results where available. 

PACS numbers: 31.15.ac, 31.15.ap, 34.20.Cf 



I. INTRODUCTION 



Studies of the Be^ ion are of interest due to its importance in a number of applications. First, 
the Be+ ion is used as an auxiliary ion to sympatheticahy cool other atomic or ionic species [1, 2] 
that cannot be directly laser cooled due to the lack of closed optical transitions. Second, the 
Be^ ion can combine with other atoms or neutral molecules to form molecular ions, such as 
BeH+ [3, 4, 5, 6], and the study of the long-range interaction between a Be"*" ion and atoms 
or molecules may open new routes for the study of state-elective chemical reactions relevant to 
astrophysics [3]. Third, investigations of Be"^ ion collisions with rare gases would be useful in the 
study of ion-atom Feshbach resonances [7] , pressure broadening of alkaline-earth-metal ions [8] , and 
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in studies of excitation spectroscopy of the collision [9, 10]. Fourth, since beryllium has a number 
of isotopes, studies of the Be+ ion could potentially be used to determine the nuclear charge radii 
of beryllium isotopes [11]. Finally, there is interest in studying the spectrum of the alkaline-earth- 
metal atom in high angular momentum Rydberg states. Experimental investigations have been 
made on a number of atoms [12, 13, 14, 15] with a view to determining the polarizabilitics of 
the singly ionized parent ion. These experiments measure the high (n, L) energy splitting and 
then use a polarization model to extract the polarizabilitics. Recent calculations [16] have shown 
the polarization model to be sensitive to nonadiabatic effects. The Be+ ion would be a useful 
candidate for a validation experiment since it should be less sensitive to adiabatic effects and its 
polarizabilitics can be calculated to very high precision. 

The above physical phenomena are influenced by the properties of the Be+ ion, and in particular 
the polarizabilitics and hypcrpolarizabilities. There have been several calculations of the Be"'" ion 
polarizabilitics reported in the literature. These include the work by Adelman and Szabo [17] using 
the Coulomb-like approximation, the calculation by Pipin and Woznicki [18] using the variation- 
perturbation approach with a combined Hylleraas-configuration interaction (CI) basis set, the 
calculation by Patil and Tang [19] using the valence electron binding energy to construct wave 
functions constrained to have the correct long-range asymptotic behavior, and finally the large 
basis full core plus configuration interaction (FCCI) calculations by Wang and collaborators [20, 
21, 22, 23, 24, 25]. However, there have been no calculations reported on the polarizabilitics and 
hypcrpolarizabilities for the Be+ ion excited states. 

In this paper, the polarizabilitics, hypcrpolarizabilities and some long range ion-atom disper- 
sion coefficients involving Be"*" ion are computed with two independent methods. First, oscillator 
strengths for many low-lying transitions are determined. Next, the polarizabilitics and hypcrpolar- 
izabilities for the 2 ^5 state and 2 state of Be"^ ion are computed variationally using expansions 
of the wave functions in Hylleraas bases. The same set of long-range parameters are also com- 
puted using a fixed core plus semi-empirical polarization potential to describe the valence electron. 
The agreement between the two different calculations will be seen to be excellent. The long-range 
dispersion interactions between Bc^ ion and the H, He or Li atoms are given and once again the 
agreement between the two sets of calculations is excellent. Furthermore, we compute all the pa- 
rameters needed to define a Be"*" polarization potential (including terms up to r~*) to describe the 
high-L Rydberg states of beryllium, where r is the ion-electron distance. All results of this paper 
are given in atomic units {e = h = irie = 1)- 
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II. THEORY AND METHOD 
A. Hylleraas variational method 

The calculations for Be^ ion are very similar in style to those for Li [26, 27]. In the center of 
mass frame, the nucleus is taken as the reference particle 0, with mass mo and charge qo, rj is the 
electron-nucleus distance, and i = 1,2,3. The nonrelativistic Hamiltonian of this system can be 
written in the form 

^0 = -E^v?-- E v..v,+,oE-+ E — ' (1) 

^ 2ni mo ^ ri 

where rij = |rj — r^l is the distance between electrons i and j, qi are the charge of the three 
electrons, and jii = miiJio / {rrii + uiq) is the reduced mass between the ith. electron and the nucleus. 
In the present paper, all the calculations are done in the infinite nuclear mass approximation. 

Significant progress has been made recently in variational calculations for three electron systems 
by using multiple basis sets in Hylleraas coordinates [28, 29]. These have the functional form, 

^ = rf rfrf r{-riVf,^e---/^^-^'^3^JM, 2, 3) , (2) 

where y^^^i^Y^i is ^ vector-coupled product of spherical harmonics to form an eigenstate of total 
angular momentum L and projection M^, and x(l,2,3) is the three-electron spin function. The 
variational wave function is a linear combination of anti-symmetrized basis functions (f). With some 
truncations to avoid potential numerical linear dependence, all terms in Eq. (2) are included such 
that 

ji + 32 + h + ii2 + i23 + hi (3) 

where O is an integer, and the convergence for the energy eigenvalue is studied by progressively 
increasing $7. 

For the He atom, taking the nucleus as the reference particle 0, the electron is labeled as particle 
1 and the other electron is labeled as particle 2. The wave functions are expanded in terms of the 
explicitly correlated basis set in Hylleraas coordinates: 

= r{T{T'l^e-''^^-^'-yfX{^^M) ■ (4) 
For the hydrogen atom, we use a basis set of form 

^ = rV/^^/^Lf+^)(/3r):y,^(r), (5) 
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where L„ is the generahzed Laguerre polynomial and the parameter /? is chosen to be 

/3 = 2/(^ + 1). This basis set has been proven to be numerically stable as the size of the basis set 
is increased. 

B. Single electron model 

The detailed description of the procedure used to construct the frozen-core Hamiltonian and 
the semi-empirical polarization potential can be found in previous works by Mitroy and collabora- 
tors [16, 30, 31]. Accordingly, only the briefest description is given here. 

Initially, a Hartree-Fock calculation of the Be"^ ion ground state was performed. The core 
orbitals were then fixed, and a semi-empirical core polarization potential was added to the Hamil- 
tonian. The core dipole polarizability was taken to be 0.0523 [30]. The cutoff parameters in the 
semi-empirical core polarization potential were po = 0.941, p\ = 0.895, p2 = 1.200 with all other 
PL set to 1.00 [30]. These values were chosen to reproduce the binding energies of the low-lying 
states. 

The low-lying states and pseudo-continuum states were obtained by diagonalizing the fixed core 
Hamiltonian in a large basis of Laguerre type orbitals. Typically, there were about 50 orbitals for 
each value of the valence angular momentum. This is large enough to eliminate the basis as a 
significant source of error. 

The multipole matrix elements and oscillator strengths were computed with a modified transi- 
tion operator [30]. The adjustable parameter in the modified multiple operator was set to p = 1.00. 
Core excitations arc included in the dispersion parameter calculation. Oscillator strength distribu- 
tions were constructed from independent estimates of the core polarizabilities (a2 = 0.01532 [30] 
and as = 0.01125 [32]). The calculations using this approach are termed the Hartree-Fock plus 
core polarization (HFCP) model. 

The dispersion coefficient calculations involving H and He used matrix element lists that were 
generated using basis functions that were similar (but not identical) in construction to Eq. (4) and 
(5). The matrix element lists for the Li atom came from a one electron model as described in 
Ref. [33]. 
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C. Polarizabilities 

Using perturbation theory, the adiabatic long-range interaction potential for an ion-atom system 
up to the term can be written [34] as 

e=l n=3 

where Qa = Qi is the total charge of the ion a, af'^ is the 2^-polc static polarizabilities for 
the atom b, R is the distance between the ion and the atom, and the C2n parameters are the 
dispersion coefficients. The first term in Eq. (6) is the polarization interaction, which does not 
lead to a frequency shift between the different states of the ion. The second term is the dispersion 
interaction, which can lead to a frequency shift between two different ion states when the ion is 
immersed in a buffer gas. Prom Eq. (6) , we see that the establishment of the interaction potential 
Vab{R) accurate to R~^ requires the static polarizabilities ai, 02, and 03 for the atom 6, and the 
dispersion coefficients Cq and Cs between ion a and atom b. Once we have obtained the oscillator 
strength spectra between the ground states and the intermediate states for the H, He, Li, and 
Be"*" ion systems, we can calculate the Be"*" ion polarizabilities, and the dispersion coefficients for 
various combinations of these particles. The detailed derivation of the formulae for calculating 
the polarizabilities, hyperpolarizabilities and dispersion coefficients in Hylleraas coordinates can 
be found in the Ref. [26]. 

For the high-L Rydberg states of an atom or an ion, where a single electron is in highly excited 
state but still moves in the field of a spherically symmetric core, the polarization interaction between 
the core and a single Rydberg electron leads to the effective potential [12, 13, 16, 35, 36], 

^eff(r) = — ^ - - ^ - ^8 + •••, (7) 

where r is the radial coordinate of the Rydberg electron and the coefficients An are properties of 
the free ion core. Comparisons with high precision variational calculations of the Rydberg states 
for the few-body systems. He [37], Li [21], and H2 [38] have been made. This functional form has 
also been used to analyze the fine structure spectrum of the Rydberg states of neutral Mg and 
Ba, resulting in estimates of the dipole polarizabilities of the alkali-earth-metal ions Mg+ and Ba"^ 
ground states [12, 14, 31]. 

According to definitions given previously [31, 35, 39], the leading coefficient A^ is half the size 
of the static dipole polarizability, 

^4 = f, (8) 
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with 

_87r^ KnoO||ri||nl)p 

The notation \nL) indicates the intermediate state with main quantum number n and angular mo- 
mentum number L. and Ti is the dipole transition operator, which satisfies the general expression 
for the 2^-pole transition operator in the center of mass frame, 

3 

Te = Y,<iyiYeoiiri). (10) 

i=l 

The next term Aq is composed of two separate terms, 

Ae = ^^, (11) 

where a2 is the quadrupole polarizability and /3i is the first order nonadiabatic correction to dipole 
polarizability. They are defined as 

47r^ KnoO||ri||nl)p 



9 ^ (^„-£^o)^ 
The term A-; also comes from two parts, namely, 

A, = -'Jl±^, (14) 

where q is the charge on the core, rj arises from third order perturbation theory and it can be 
expressed as 

327r\/l07f (noO||ri||ml)(ml||ri||n2)(n2||r2||noO) 
" 225 {E^-Eo){En-Eo) 



167rV67r ^-^ (noO||Ti ||ml)(ml||T2 ||nl)(nl||ri ||noO) 

135 £^ (S^-i?o)(i?„-i^o) ' ^ ^ 

and (5i is the second order nonadiabatic correction to the dipole polarization energy, 

47r^ KnoO||ri||nl)p 

Quite a few terms contribute to As, 

as - 15/32 + e - ai/^i + 72(5i 
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The octupole polarizability 0:3 is computed by 

«3=49i: , (18) 

while f32 comes from the first nonadiabatic correction part to the quadrupole polarization energy, 
and is defined as 

47r^ |(no0||r2||n2)p 



n 



The term e is defined 



^ _ 327r2 ^ |(noO||ri||ml)(ml||Ti||nO)|2 

m,n 



81 (ii;™-i5o)2(i^n-i5o) 



647r2^ |(noO||ri||ml)(ml||ri||n2)p 

405 ^ (£;^-Eo)2(^n-^o) ■ ^ ' 



m.n 



The hyperpolarizability 70 of the ground state, and coefficients e, ai, /3i are related by the identity 
70 = 12(e — ai/3i). The last term A^l is nonadiabatic in origin and defined by 

^. = ^. (21) 

For excited states, the working expressions for the evaluation of the polarizabilitics ai, , a2, 
and Q!3 and hyperpolarizabilities 70 and 72 are given in Tang et al [26]. These expressions are quite 
lengthy and not reproduced here. 



D. Dispersion interactions 

The dispersion interaction, between two atoms, A and B, for the Hylleraas wave functions was 
calculated from tabulated lists of matrix elements using sum rules [26] . For the case of Cq involving 
two S-state atoms, one can write the sum rule 



The sum is over all states of P°-symmetry. The absorption oscillator strength, ,for a dipole 
transition from — n, with an energy difference of AEmO = Eq — Em, is most conveniently defined 
in the present context [27, 30] as 

= 3(2Lo + 1) • 
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In this expression for the HFCP method, Lq is the ground state orbital angular momentum, and 
rC^(f) = y'^rY^(f) is the operator for a dipole transition just for one electron. Explicit expres- 
sions for Cs, and Cio, and for other symmetries can be found in [26]. 

The dispersion interaction calculations for the HFCP wave functions utilized the completely 
general procedures outlined by Zhang and Mitroy [40, 41]. These calculations utilize sum rules 
involving lists of reduced matrix elements of the multipolc operator r^C^ multiplied by angular 
recoupling factors. They are effectively equivalent to Eq. (22) and the expressions in [26] despite 
the differences how the calculations are actually carried out. 



III. RESULTS AND DISCUSSION 
A. The energies and oscillator strengths of Be"*" ion 

Table I shows the convergence study for the nonrelativistic energy of the Be"*" ion ground state 
as the size of the Hylleraas basis set is enlarged. The ratio R{^) is defined by 

^ ' E{n)-E{n-i) ^ ' 

The extrapolation was done by assuming that the ratio would stay constant around R{Vt) = 4.839. 
It is clear from Table I that the energy converges to high precision as the number of terms is 
progressively increased. The final nonrelativistic energy is accurate to about 11 significant figures 
for the ground state. The uncertainty in the final energy is set to be equal to the extrapolation 
correction. The nonrelativistic theoretical energies for other low- lying states of Be"*" ion are tab- 
ulated and compared with experimental binding energies in Table II. Binding energies from the 
HFCP calculations are also listed. The Hylleraas binding energies are generally in good agreement 
with the experimental binding energies. There are two effects which change the binding energies, 
finite mass and relativistic effect. The finite mass effect will be about 0.001% and will probably 
act to decrease the magnitude of the binding energy. The largest discrepancy for the 2 state is 
due entirely to relativistic effects (see Table III of [29]), which contributes about 0.007% of the 
final value. 

Table III lists the oscillator strengths for a number of the Be"*" ion dipole transitions involving 
low-lying states. The final values for the Hylleraas calculations are obtained with an extrapolation 
procedure similar to that for the energy. Once again the uncertainty in the oscillator strength is 
assigned to be equal to the magnitude of the extrapolation correction. It is evident from Table HI 
that there is excellent agreement between the Hylleraas and HFCP calculations. The largest 
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discrepancy between the two calculations is only about 0.1% (for the 3^P ^ 3^D transition). The 
oscillator strengths from the Hylleraas calculation could be used to improve the National Institute 
of Standards and Technology (NIST) tabulations [42, 43]. 

The Hylleraas oscillator strength for the 2^S ^ transition is accurate to about 7 significant 
figures and is compatible with an earlier Hylleraas calculation by Yan et al. [44]. Some earlier large 
scale ab-initio calculations also gave oscillator strengths that are compatible with the present Hyller- 
aas calculation for this transition. These include the multi-configuration Hartree-Fock (MCHF) 
calculation of Godefroid et al. [45] and the FCCI calculation of Chung et al [22] . While the present 
oscillator strengths are reported with 7 significant digits, finite-mass and relativistic effects that 
are not included in the calculation could conceivably alter the oscillator strengths beginning at the 
fifth digit. 

There have also been some high precision oscillator strengths reported for the 2^P — > 3^D 
transition. The MCHF value of Godefroid et al. [45] and the FCCI value of Qu et al. [23] agree 
with the Hylleraas calculation to better than four digits. Qu et al. [24, 25] have also reported 
GI calculations for the 2^5 — >■ 3^P and 3^1? — > 4^F transition and again there is agreement to 
better than four significant digits. The experimental oscillator strengths listed in the Table were 
measured using beam-foil spectroscopy [46] . These oscillator strengths have only low precision and 
cannot discriminate between the higher quality theoretical estimates. 

B. Polarizabilities of the Be+ ion 

The convergence properties of the static dipole polarizability ai, and hyperpolarizability 70 for 
the Bc+ ion ground state are presented in Table IV. Both of them have converged to five significant 
figures. The extrapolation was done by assuming that the ratio between two successive diS'erences 
would stay constant as the basis size increased towards infinity. The uncertainty in the final value 
is set equal to the magnitude of the extrapolation correction from the explicitly calculated value 
computed with the basis of largest dimension. 

Table V shows the convergence as a function of basis size for the scalar and tensor dipole po- 
larizabilities a\, aj, and the hyperpolarizabilities 70, 72 for the first excited state of Be+ ion. The 
intermediate sums in this case have contributions from doubly excited unnatural parity states with 

= and = 2° (e.g. the unnatural parity 1*^ state has two i = 1 electrons coupled to a total 
angular momentum of L = 1). The contribution from the unnatural parity states is usually small. 
For example, the unnatural parity states contribution of 0.020616 to ai is about 1%. The scalar 
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dipole polarizability is converged to five significant digits, and the tensor dipole polarizability aj is 
converged to six digits. The hyperpolarizabihties, 70 and 72 are accurate to the six and five signif- 
icant figures respectively. There was no major numerical cancellation in the hyperpolarizabihties 
in our calculation, in contrast to the situation that prevails for the Li hyperpolarizability, which 
suffers severely from cancellations in the different parts of the calculation [26] . 

There have been a number of accurate calculations of the multipolar polarizabilities for the Li 
atom in its ground or lowest excited states [26, 33, 50]. However, there have been fewer polariz- 
ability calculations for the Be+ ion in its ground and lowest energy excited state. Table VI gives a 
comparison between the present results and previous calculations for the static polarizabilities of 
Be+ ion in the 2 and 2 states. One of the most notable features of the Table is the very good 
agreement between the Hylleraas and HFCP multipole polarizabilities. The overall level of agree- 
ment is at the 0.1% level. The one exception was the static dipole polarizability of the 2^P state, 
but it should be noted that this polarizability is small due to cancellation between different terms 
in the oscillator strength sum. The polarizability of Be"'"(2^i-*) state is relatively small because the 
Be"*" (2 — s- 2 ^5) oscillator strength is negative while all the other Be"*" (2 — > n ^L) /-values are 
positive. The net effect of the cancellations is a reduction in the overall size of the polarizability 
by a factor of five. 

The Coulomb approximation polarizability [17] and the asymptotically correct wave function 
polarizabilities [19] achieve about 1-2% accuracy in ai and a2- This is noticeably worse than any 
of the other polarizabilities listed in the Table. The older Hylleraas-typc calculations by Pipin and 
Woznicki [18] used the variation-perturbation approach to estimate the polarizabilities (as opposed 
to oscillator strength sum rules). Their value of ai = 24.5 [18] is compatible with the present value 
but not nearly as precise as the Hylleraas polarizability. The variation-perturbation calculations 
using the FCCI wave function [20, 21] gave polarizabilities that agree with the present Hylleraas 
polarizabilities to better than 0.01%. 

The present results for the Be+(2^P) state are the only results reported for the higher order 
polarizabilities and the hyperpolarizabihties. Previously Merawa and Rerat reported the calcu- 
lations for ai and of for the Be"'"(2^P) state using the time-dependent gavigc-invariant method 
(TDGI) [51], but the underlying structure model for this approach is less accurate than the present 
calculations and we do not include their numerical values in Table VL The overall level of agree- 
ment between the Hylleraas and HFCP calculations is very impressive when it is considered that 
there are significant numerical cancellations in the calculation of ai that lead to a small value. 



11 



C. Effective potential for beryllium Rydberg state 

Recently there have been a number of investigations of ion polarizabihties based on the inter- 
pretation of resonant excitation stark ionization spectroscopy (RESIS) [12, 14, 15, 16]. The energy 
sphtting of adjacent Rydberg levels with AL = 1 is used to determine the parent ion polarizabih- 
ties. One recent finding has been an increased appreciation of the importance of nonadiabatic and 
higher order polarizability terms proportional to and in the interpretation of the RESIS 
spectra. 

Tabic VII summarizes all the parameters necessary to define the polarization series given by 
Eq. (7) for the Be"*" ion ground state. The data are presented since the neutral beryllium series 
represents an ideal system upon which to validate the underlying assumptions used in the analysis of 
the RESIS experiment. The nonadiabatic effects are strong [the nonadiabatic dipole polarizability 
of /3i = —81.78175(1) dominates the quadrupole polarizability of a2 = 53.7659(2) in the evaluation 
of Aq = —218.4622(1)] and all the "Hylleraas" polarizabihties listed in Table VII would have an 
overall level of precision better than 0.1%. Although no experiment has been done, a RESIS 
experiment upon neutral beryllium would provide a stringent test on the ability of an analysis 
based on Eq. (7) to extract polarizabihties from a typical RESIS spectrum. 

D. The long-range dispersion coefficients 

Table VIII and Table IX list the long-range dispersion coefficients for the Be+ ion interacting 
with the H, He, and Li atoms. Table VIII lists dispersion coefficients when both atoms or ions are 
in their ground states. Table IX gives dispersion coefficients when one of the systems in an excited 
state. All of the dispersion coefficients have been calculated independently using the Hylleraas and 
HFCP wave functions. Besides the wave functions, the procedures used to combine the lists of 
matrix elements were completely independent. As far as we know, the data listed in Tables VIII 
and Tables IX are the only dispersion coefficients published for these systems. 

The level of agreement between the two sets of C„ values is generally excellent. For example, 
the largest difference between any of the dispersion constants listed in Table VIII is only 0.06%, 
occurring for the Be+{2'^S)-Li{2'^S) value of Cg. 

The high level of agreement also occurs for the Cn values listed in Table IX, the only case of a 
greater than 1% difference occurring for the Be+(2 ^P)-Li(2 ^5) dimer. In this case, the roughly 3% 
disagreement occurs as a result of the previously mentioned cancellations in the oscillator strength 
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sum for the Be+ (2 ^P) polarizability. The net effect of the cancellations is a reduction in the overall 
size of the dispersion constants by a factor of about 100. For example, the first term (A = 0) of 
Eq. (52) of Tang et al. [26] was -115.545 while the second term( A = 2) was 117.3969. 

We do not list dispersion coefficients for the state combinations that allow Penning or associative 
ionization (this occurs when the excitation energy of one atom is sufficient to cause ionization in 
the other atom). When this is possible, there is a singularity in the energy denominator of the 
oscillator strength sum rules which makes it problematic to achieve convergence. 

All the values in Table VIII and Table IX provide an important benchmark for the accurate 
determination of the interaction potentials between Be"^ ion and the H, He, or Li atoms. The BeH+ 
ion is one of the few molecular ions [3] for which a potential curve could be computed with an 
explicitly correlated wave function. Hence the present values for the Be+-H long-range interaction 
could be used to help construct a very accurate global potential surface for this system. 

IV. CONCLUSION 

Fully correlated Hylleraas variational wave functions have been used to determine definitive 
values for the oscillator strengths, polarizabilities and hyperpolarizabilities for the Be"*" ion 
ground state and the 2 excited state. The Hylleraas results for the polarizabilities of the 2 "^S 
state improve the accuracy of previous values by more than one order of magnitude. Complemen- 
tary calculations using a semi-empirical method have also been done. The high level of agreement 
between the two calculations at the 0.1% level of precision attests to the utility of carefully formu- 
lated effective potential approaches, which can give good descriptions of atomic structure with low 
computational expense. 

The long-range dispersion coefficients for the Be"*" ion interacting with a H, a He, or a Li atom 
have been evaluated. The polarizabilities and dispersion coefficients provide reliable references for 
the description of ion-atom collisions involving Be''" ion and also for high precision calculations of 
the potential curves between the Be+ ion and atoms such as H, He or Li. 

In addition, all the parameters of the effective polarization potential for the Be"*" ion up to the 
term have been obtained. These parameters are extremely useful in the description of high-L 
Rydberg states of beryllium, and could be used in future experiments to determine the ionization 
potential of beryllium and also to describe the fine structure of beryllium atom Rydberg series. In 
addition, the present calculations lay the foundation for the further investigations of relativistic 
and QED effects on the polarizabilities and other properties of the Be+ ion. 



13 



Acknowledgments 

This work was supported by NNSF of China under Grant No. 10674154 and by the National 
Basic Research Program of China under Grant No. 2005CB724508. Z.-C.Y. was supported by 
NSERC of Canada and by the computing facihties of ACEnet, SHARCnet, WestGrid, and in 
Part by the CAS/SAFEA International Partnership Program for Creative Research Teams. J.F.B 

was supported in part by the U.S. NSF through a grant for the Institute of Theoretical Atomic, 
Molecular and Optical Physics at Harvard University and Smithsonian Astrophysical Observatory. 
J.M. and J.Y.Z would like to thank the Wuhan Institute of Physics and Mathematics for its 
hospitality during their visits. 



T. Rosenband et al, Science. 319, 1808 (2008). 

T. Rosenband et al, Phys. Rev. Lett. 98, 220801 (2007). 

B. Roth, P. Blythe, H. Wenz, H. Daerr and S. Schiller, Phys. Rev. A 73, 042712 (2006). 

B. Roth, P. Blythe, H. Daerr, L. Patacchini and S. Schiller, J. Phys. B 39, S1241 (2006). 

J. B. Roos, M. Larsson, A. Larson and A. E. Orel, Phys. Rev. A 80, 012501 (2009). 

L. F. Errea, B. Herrero, L. Mendez, 1. Rabadan and P. Sanchez, J. Phys. B 27, L753 (1994). 

Z. Idziaszek, T. Calarco, P. S. Julienne and A. Simoni, Phys. Rev. A 79, 010702(R) (2009). 

P. S. Barklem and B. J. O'Mara, Mon. Not. R. Astron. Soc. 311, 535 (2000). 

J. 0. Olsen and T. Andersen, J. Phys. B 8, L421 (1975). 

N. Andersen, T. Andersen, J. 0. Olsen and E. H. Pedersen, J. Phys. B 13, 2421 (1980). 
W. Nortershauser et al, Phys. Rev. Lett. 102, 062503 (2009). 

E. L. Snow, M. A. Gearba, R. A. Komara, S. R. Lundeen and W. G. Sturrus, Phys. Rev. A 71, 022510 
(2005). 

E. L. Snow and S. R. Lundeen, Phys. Rev. A 75, 062512 (2007). 
E. L. Snow and S. R. Lundeen, Phys. Rev. A 77, 052501 (2008). 

B. J. Lyons and T. F. Gallagher, Phys. Rev. A 57, 2426 (1998). 
J. Mitroy, Phys. Rev. A 78, 052515 (2008). 

S. A. Adelman and A. Szabo, J. Chem. Phys. 58, 687 (1973). 
J. Pipin and W. Woznicki, Chem. Phys. Lett. 95, 392 (1983). 
S. H. Patil and K. T. Tang, J. Chem. Phys. 106, 2298 (1997). 
Z. W. Wang and K. T. Chung, J. Phys. B 27, 855 (1994). 

C. Chen and Z. W. Wang, J. Chem. Phys. 121, 4171 (2004). 

K. T. Chung, in Proceedings of the Vlth International Conference on the Physics of Highly Charged 
Ions, edited by P. Richard, M. Stockl, C. L. Cocke and C. D. Lin (Manhattan, Kansas, 1992), AIP 



14 



Conf. Proc. 274 (AIP, New York, 1993). Z.-W. Wang, X.-W. Zhu and K. T. Chung, Phys. Rev. A 46, 
6914 (1992). 

[23] L. H. Qu, Z. W. Wang and B. W. Li, J. Phys. B 31, 3601 (1998). 

[24] L. H. Qu, Z. W. Wang and B. W. Li, Opt. Comm. 162, 223 (1999). 

[25] L. H. Qu, Z. W. Wang and B. W. Li, Eur. Phys. J. D. 5, 175 (1999). 

[26] L. Y. Tang, Z. C. Yan, T. Y. Shi and J. F. Babb, Phys. Rev. A 79, 062712 (2009). 

[27] Z. C. Yan, J. F. Babb, A. Dalgarno and G. W. F. Drake, Phys. Rev. A 54, 2824 (1996). 

[28] M. Puchalski and K. Pachucki, Phys. Rev. A 78, 052511 (2008). 

[29] Z. C. Yan, W. Nortershauser and G. W. F. Drake, Phys. Rev. Lett. 100, 243002 (2008). 

[30] J. Mitroy and M. W. J. Bromley, Phys. Rev. A 68, 052714 (2003). 

[31] J. Mitroy and M. S. Safronova, Phys. Rev. A 79, 012513 (2009). 

[32] R. P. McEachran, A. D. Stauffer and S. Greita, J. Phys. B. 12, 3119 (1979). 

[33] J. Y. Zhang, J. Mitroy and M. W. J. Bromley, Phys. Rev. A 75, 042509 (2007). 

[34] J. Y. Zhang and Z. C. Yan, J. Phys. B: At. Mol. Opt. Phys. 37, 723 (2004). 

[35] R. J. Drachman, Phys. Rev. A 26, 1228 (1982). 

[36] R. J. Drachman and A. K. Bhatia, Phys. Rev. A 51, 2926 (1995). 

[37] Z. C. Yan, Phys. Rev. A 62, 052502 (2000). 

[38] P. L. Jacobson, D. S. Fisher, C. W. Fehrenbach, W. G. Sturrus and S. R. Lundeen, Phys. Rev. A 56, 
R4361 (1997). 

[39] G. W. F. Drake, Adv. At. Mol. Opt. Phys. 31, 1 (1993). 

[40] J. Y. Zhang and J. Mitroy, Phys. Rev. A 76, 022705 (2007). 

[41] J. Mitroy and J. Y. Zhang, Phys. Rev. A 76, 032706 (2007). 

[42] G. A. Martin and W. L. Wiese, Phys. Rev. A 13, 699 (1976). 

[43] W. L. Wiese, M. W. Smith and B. M. Glennon, NSRDA-NBS 4, (1966). 

[44] Z. C. Yan, M. Tambasco and G. W. F. Drake, Phys. Rev. A 57, 1652 (1998). 

[45] M. Godefroid, C. F. Fischer and P. Jonsson, J. Phys. B 34, 1079 (2001). 

[46] T. Andersen, K. A. Jessen and G. S^rensen, Phys. Rev. 188, 76 (1969). 

[47] Y. Ralchenko, A. E. Kramida, J. Reader and NIST ASD Team (2008). NIST Atomic Spectra Database, 
Version 3.1.5. URL http://physics.nist.gov/asd3. National Institute of Standards and Technology, 
Gaithersburg, MD. 

[48] S. Bashkin and J. O. Stoner, Jr., Atomic energy levels and Grotrian diagrams (American Elsevier, New 
York, 1975). 

[49] G. W. F. Drake, Handbook of Atomic, Molecular, and Optical Physics (American Institute of Physics, 
New York, 1996). 

[50] W. R. Johnson, U. I. Safronova, A. Derevianko and M. S. Safronova, Phys. Rev. A 77, 022510 (2008). 
[51] M. Merawa and M. Rerat, J. Chem. Phys. 108, 7060 (1998). 



15 



TABLE I: Convergence of the Hylleraas calculation of the nonrel- 
ativistic energy (in atomic units) for the 2 ^5 state of Be+ ion. 



n 


No. of t(^rms 


E{0) 


E{0)-EiO-r) 


R(fi) 


8 


1589 


-14.324 763 166 358 


-0.000 000 051 721 


8.211 


9 


2625 


-14.324 763 174 596 


-0.000 000 008 238 


6.278 


10 


4172 


-14.324 763 176 309 


-0.000 000 001 713 


4.809 


11 


6412 


-14.324 763 176 663 


-0.000 000 000 354 


4.839 


ExtrM.p. 




14.;-!24 Ua 7.S()(7:-!) 







TABLE II: Theoretical nonrelativistic energies and experimental 
energies of the low- lying states for the Be+ ion, (in atomic units). 
The numerical uncertainty in the theoretical energies are given in 
brackets. The experimental valence binding energies are taken from 
the National Institute of Standards database [47, 48] . The ground- 
state energy for the Be^+ ion was taken from Ref. [49]. 



State 


Hylleraas 


SHyUcraas-S(Be2+) 


HFCP 


Experiment 


Ref. [29] 


2^S 


-14.324763176736(73) 


-0.669196938312 


-0.669250 


0.669247 


0.66924793(2) 


2 2p 


-14.17933329329(24) 


-0.52376705486 


-0.523755 


0.523769 


-0.52376988(2) 




-13.9227892683(5) 


-0.2672230298 


-0.267189 


-0.267233 


-0.26723367(3) 


3 2p 


-13.8851502898(5) 


-0.2295840513 


-0.229527 


-0.229582 




3^D 


-13.87805405934(36) 


0.22248782091 


-0.222482 


-0.222478 






-13.7987166133(8) 


0.1431503748 


0.143131 


0.143152 






-13.780581705614(80) 


0.125015467190 


0.125015 


0.125008 




5^ 


-13.735568352173(16) 


-0.080002113749 


-0.080002 


-0.079997 





TABLE IIL Dipole oscillator strengths for the selected transitions 
of Be"*" ion. 



Transition Hylleraas 



HFCP NIST [42, 43] Exp. [46] 



Other theory 



2^S^2^P 0.49806736(6) 0.4985 0.505 

2 25^32p 0.08316525(18) 0.0828 0.0804 

2 2p^325 0.06434157(29) 0.0643 0.0665 

2 2p^425 0.01021583(30) 0.0102 0.010 



0.54(3) 



0.048(5) 



0.498067381(25) Hylleraas [44] 
0.49813 FCCI [22]; 0.49807 MCHF [45] 
0.08136 FCCI [25] 
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2^P^3^D 0.6319828(11) 0.6321 0.652 

S'^S^S^P 0.8297696(15) 0.8307 

3 2p^3 2£, 0.08103350(17) 0.0804 

32p^425 0.1345245(13) 0.1346 

32^)^42^ 1.01460194(11) 1.0146 1.01 

A'^F — r>"-^G :i.:-')4r).S712(i(12) l.:Ur,:] 



0.63199 MCHF [45] 
0.63197 FCCI [23] 



0.66(3) 1.0146 FCCI [24] 



TABLE IV: Convergence of the dipole polarizability ai and hyper- 
polarizability 70 for the Be^ ion ground state (in atomic units). 
The number of intermediate states of a given angular momentum 
are denoted as Ns, Np and Nn- 



iNs,Np) 


ai 


{Ns,Np,Nd) 


70 


(1589,1174) 


24.495 332 


(1589,1174,1174) 


115 21.13 20 


(2625,2091) 


24.496 067 


(2625,2091,2091) 


115 21.31 84 


(4172,3543) 


24.496 408 


(4172,3543,3543) 


-115 21.31 96 


(6412,5761) 


24.496 522 


(6412,5761,5761) 


-115 21.27 68 


Extrap. 


24.496 6(1) 


Extrap. 


-115 21.30(3) 



TABLE V: Convergence of the dipole polarizabilities ai, of and 
the hyperpolarizabilities 70, 72 for the Be+ ion in the 2 state (in 
atomic units). The number of natural parity intermediate states 
of a given angular momentum are denoted as Ns, Np, Nn and 
Np. The number of unnatural parity intermediate states of a given 
angular momentum are denoted as Npi and Nn'- 



{Ns,Np,Nn,Np,) 


ai 




{Ns,Np,Nn,Np,Np,,ND') 


70 


72 


(1589,1174,1174,1106) 


2.024 6197 


5.856 054 26 


(1589,1174,1174,1248,1106,1428) 


109 11.665 61 


-7372.0 881 


(2625,2091,2091,2002) 


2.024 7235 


5.856 019 68 


(2625,2091,2091,2307,2002,2640) 


109 13.221 87 


-7373.4 188 


(4172,3543,3543,3413) 


2.024 7465 


5.856 014 25 


(4172,3543,3543,4051,3413,4587) 


109 13.576 50 


-7373.6 683 


(6412,5761,5761,3413) 


2.024 7537 


5.856 013 46 


(6412,5761,5761,6806,3413,4587) 


109 13.572 18 


-7373.5 698 


Extrap. 


2.024 76(1) 


5.856 012(1) 


Extrap. 


109 13.57(1) 


-7373.6 1(5) 
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TABLE VI: Comparisons of the static polarizabilities and hyper- 
polarizabilities (in atomic units) for the 2 and 2 states of Be+ 
ion. 











2 25 state 






Method 




ai 






70 




Coulomb approximation [17] 


24.77 










Variation-perturbation 


Hylleraas CI [18] 24.5 










Asymptotic correct wave function [19] 


24.91 


53.01 


400. / 






Variation-perturbation 


FCCI [20, 21] 


24.495 


53.774(24) 


465.79(11) 






HFCP 




24.493 


53.760 


465.77 


-11511 




Hylleraas 




24.4966(1) 53.7659(2) 


465.7621(1) 


-11521.30(3) 












2 2p state 






Method 




Oil 




0:2 


as 


70 72 


HFCP 




2.028 


5.835 


62.313 


1208.8 


10996 -7450.4 


H\-l](^ra,;is 




2.0247(i(l) 


r).8.",(i012('l^ 


1 (i2. 2840(1) 


1207.812(2) 


1091:-!.. ")7(1) 7.'-!7.'-!.()l(r)) 



TABLE VII: The polarizability parameters and coefficients An of 
the polarization potential Eq. (7). 



Method 


ai 


Pi 


Q!2 


/32 


Q!3 


Si 


V 


e 


Hylleraas 


24.4966(1) 


81.78175(1) 


53.7659(2) 


58.1169(1) 


465.7621(1) 


279.16401(2) 


917.569(1) 


1043.27(1) 


HFCP 


24.494 


81.751 


53.760 


58.106 


465.77 


278.94 


917.37 


1038.7 




A4 




A7 


^8 


Arl 








Hylleraas 


12.24830(5) 


-218.4622(1) 


905.445(2) 


9366.857(1) 


1004.99051(8) 








HFCP 


12.247 


-218.37 


-904.99 


9357.2 


1004.2 









TABLE VIII: The long-range dispersion coefficients Ce and Cg for 
the Be+ ion ground state interacting with the H, He, and Li atom. 
The first row for each system came from Hylleraas wave functions 
while the second row were computed with HFCP wave functions. 



System Cq Cg 

Be+{2'^S)-ll{l^S) 18.8314(1) 371.675(5) 

18.829 371.62 
Be+(2 2S')-He(li5) 6.9811(1) 120.425(3) 



6.979 120.44 
Be+(2 25)-He(2i5) 621.577(2) 41371.9(1) 

621.52 41370 
Be+(2 25)-He(2 3S') 400.289(3) 19753.9(2) 

400.26 19753 
Be+(2 2S')-Li(2 25) 286.75(1) 11991.1(2) 

286.82 11998 



TABLE IX: The long-range dispersion coefficients Ce, and Cg for 
the Be"*" ion interacting with H, He, and Li atom for the atomic 
states with a combined angular momentum of L = 1. The Mb 
column denotes the total magnetic quantum number of the system. 
The first row for each system came from Hylleraas wave functions 
while the second row were computed with HFCP wave functions. 



System 


Mb Ce 




Be+(2 2p)-H(l2S') 





38.53656(1) 


1291.81(1) 







38.571 


1292.7 


Be+(2 2p)-H(l25) 


±1 


17.12124(1) 


172.578(1) 




±1 


17.128 


172.67 


Be+(2 2p)-He(li5) 





11.64784(1) 


438.3026(6) 







11.660 


438.63 


Be+(2 2p)-He(li5) 


±1 


6.028208(7) 


34.0606(3) 




±1 


6.030 


34.113 


Be+(2 25)-Hc(2'^P) 





791.311(4) 


93960.4(9) 







791.52 


93960 


Be+{2^S)-Re{2^P) 


±1 


373.997(3) 


3547.53(8) 




±1 


374.25 


3553.1 


Be+(2 2p)-Li(2 2S') 





-326.06(8) 


-3.613643(6) xlO^ 







-325.74 


-3.52050x10^ 


Be+(2 2p)-Li(2 25) 


±1 


2.32(4) 


-1.214626(3) xlO^ 




±1 


2.3964 


-1.1836x10^ 


Be+(2 2S')-Li(2 2p) 





925.279(2) 


1.043315(3) xlO^ 







925.49 


1.0438 xlO^ 


Be+(2 2S')-Li(2 2p) 


±1 


420.470(1) 


3804.79(4) 




±1 


420.53 


3806.7 



